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Compton scattering of a laser beam with a relativistic electron beam has been used to generate 
intense, highly polarized and nearly monoenergetic x-ray or gamma-ray beams at many facilities. 
The ability to predict the spatial, spectral and temporal characteristics of a Compton gamma-ray 
beam is crucial for the optimization of the operation of a Compton light source as well as for the 
applications utilizing the Compton beam. In this paper, we present two approaches, one based upon 
analytical calculations and the other based upon Monte Carlo simulations, to study the Compton 
scattering process for various electron and laser beam parameters as well as different gamma-beam 
collimation conditions. These approaches have been successfully applied to characterize Compton 
gamma-ray beams, after being benchmarked against experimental results at the High Intensity 
Gamma-ray Source (HI7S) facility at Duke University. 
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I. INTRODUCTION 

Compton scattering of a laser beam with a relativistic 
electron beam has been successfully used to generate in- 
tense, highly polarized and nearly monoenergetic x-ray 
or gamma-ray beams with a tunable energy at many fa- 
cilities [H-Q • These unique Compton photon beams have 
been used in a wide range of basic and application re- 
search fields from nuclear physics to astrophysics, from 
medical research to homeland security and industrial ap- 
plications 

The ability to predict the spectral, spatial and tem- 
poral characteristics of a Compton gamma-ray beam is 
crucial for the optimization of the gamma-ray beam pro- 
duction as well as for research applications utilizing the 
beam. While the theory of particle-particle (or electron- 
photon) Compton scattering, which is equivalent to the 
scattering between a monoenergetic electron beam and a 
monoenergetic laser beam with zero transverse sizes, is 
well documented in literature , there remains a need 
to fully understand the characteristics of the gamma-ray 
beam produced by Compton scattering of a laser beam 
and an electron beam with specific spatial and energy 
distributions, i.e., the beam-beam scattering. 

Study of beam-beam Compton scattering has been re- 
cently reported in 0, H|. However, the algorithms used 
in these works are based upon the Thomson scattering 
cross section, i.e., an elastic scattering of electromag- 
netic radiation by a charged particle without the recoil 
effect. For scattering of a high energy electron beam and 
a laser beam, the recoil of the electron must be taken 
into account. The Compton scattering cross section has 
been used to study characteristics of Compton gamma- 
ray beams by Duke scientists in 1990's 0, However, 



the effects of incoming beam parameters, and the effects 
of gamma-beam collimation were not fully taken into ac- 
count. 

In this paper, we present two different methods, a 
semi- analytical calculation and a Monte Carlo simula- 
tion, to study the Compton scattering process of a po- 
larized (or unpolarized) laser beam with an unpolarized 
electron beam in the linear Compton scattering regime. 
Using these two methods, we are able to characterize a 
Compton gamma-ray beam with various laser and elec- 
tron beam parameters, arbitrary collision angles, and dif- 
ferent gamma-beam collimation conditions. 

This paper is organized as follows. In Section II, we 
first review the calculation of the Compton scattered 
photon energy for an arbitrary collision angle, and then 
introduce the scattering cross section in a Lorentz in- 
variant form. Based upon this cross section, the spatial 
and spectral distributions as well as the polarization of a 
Compton gamma-ray beam are investigated in particle- 
particle scattering cases. In Section III, we discuss the 
beam-beam Compton scattering by considering effects 
of the incoming beam parameters as well as the effect 
of the gamma-ray beam collimation. Two methods, a 
semi- analytical calculation and a Monte Carlo simula- 
tion, are then presented. Based upon the algorithms of 
these methods, two computing codes, a numerical inte- 
gration code and a Monte Carlo simulation code, have 
been developed at Duke University. The benchmarking 
results and applications of these two codes are presented 
in Section IV. The summary is given in Section V. 



II. PARTICLE-PARTICLE SCATTERING 



A. Scattered photon energy 
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A review of the calculation of scattered photon energies 
in the particle-particle scattering case is in order. Fig- 
ure [T] shows the geometry of Compton scattering of an 
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electron 
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FIG. 1: Geometry of Compton scattering of an electron and a 
photon in a lab frame coordinate system (x e , y e ,z e ) in which 
the electron is incident along the z e -direction. The incident 
photon is propagating along the direction given by the polar 
angle 8i and azimuthal angle <f>i. The collision occurs at the 
origin of the coordinate system. After the scattering, the scat- 
tered photon propagates in the direction given by the polar 
angle Of and azimuthal angle <j>f. 9 P is the angle between the 
momenta of incident and scattered photons, k and k'. The 
electron after scattering is not shown in the figure. 



electron and a photon in a laboratory frame coordinate 
system (x e , y ei z e ) in which the incident electron with 
a momentum p is moving along the z e -direction. The 
incident photon with a momentum hk (h is the Planck 
constant) is propagated along the direction with angles 
(6i,<j>i). The collision occurs at the origin of the coor- 
dinate system. After the collision, the photon with a 
momentum hk' is scattered into the direction of (Qf,(j)f). 

According to the conservation of the 4- momenta before 
and after scattering, we can have 



P- 



P 



(1) 



where p = (E e /c,p) and k = (E p / c,hk) are the 4- 
momenta of the electron and photon before the scatter- 
ing, respectively; p' = (E' e /c,p?) and kl — (E g /c,hk') are 
their 4-momenta after the scattering; E e and E p are the 
energies of the electron and photon before the scattering; 
E' e and E g are their energies after the scattering; and c 
is the speed of light. Squaring both sides of Eq. (|T|) and 
following some simple manipulations, we can obtain the 
scattered photon energy as follows, 



E„ 



(1 - /3 cos 6i)E p 



(1 -/3cos0/) + (1 -cos6 p )E p /E e 



(2) 



where (3 — v/c is the speed of the incident electron rela- 
tive to the speed of light, and 8 P is the angle between the 
momenta of the incident and scattered photons (Fig. [I}. 



For a head-on collision, Oi = tt and 8 p 
can be simplified to 



E g = 



(1 + 0)E P 



(1 -j3cos6 f ) + (1 + cos9 f )E p /E e 



Of, Eq. © 
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FIG. 2: The relation between the scattered photon energy 
(in MeV) and scattering angle in an observation plane, which 
is 60 meters downstream from the collision point. The scat- 
tered photons are produced by 800 nm photons scattering 
with 500 MeV electrons. Each concentric circle is an equi- 
energy contour curve of the energy distribution of scattered 
photons. 



Clearly, given the energies of the incident electron and 
photon, E e and E p , the scattered photon energy E g only 
depends on the scattering angle Of, independent of the 
azimuth angle 4>f. The relation between the scattered 
photon energy E g and scattering angle Of is demon- 
strated in Fig. [2j In this figure, the scattered photon 
energies E g are indicated by the quantities associated 
with the concentric circles in the observation plane, and 
the scattering angles Of are represented by the radii R 
of the circles, i.e, Of = R/L, where L = 60 meters is the 
distance between the collision point and the observation 
plane. We can see that the scattered photons with higher 
energies are concentrated around the center (Of = 0), 
while lower energy photons are distributed away from the 
center. Such a relation, in principle, allows the formation 
of a scattered photon beam with a small energy-spread 
using a simple geometrical collimation technique. 

For a small scattering angle (Of <C 1) and an ultra- 
relativistic electron (7 ^> 1), Eq. ^ can be simplified 
to 



E„ 



±1 2 E P 



1 + 7 2 0f + ±1 2 E P /E e 



(4) 



where 7 = E e / (mc 2 ) is the Lorentz factor of the electron 
and mc 2 is its rest energy. When the photon is scattered 
into the backward direction of the incident photon (i.e., 
Of = 0, sometimes called backscattering) , the scattered 
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TABLE I: Relative uncertainty of the scattered photon en- 
ergy AE g /E g due to the uncertainties of various variables in 
Eq. under assumptions of ~ ir and Of ~ 0. 



Variables 


Contributions 


Approximated contributions 


E e 
E p 


Z U l+i-f 2 E p /E e ) Ee 
1 AS P 
l+i^E„/E e E v 


9 A£ t 

Be 

A Bp 
Bp 


Of 
0i 


l+i^E p /E c Af V 

-|A0? 





photon energy will reach the maximum value given by 

4 7 2 £ p 



rpmax 
E 9 



1 + ij 2 E p /E e 



(5) 



Neglecting the recoil effect, i.e., A^ 2 E p jE e < 1, Eq. j5]) 
can be reduced to the result given by the relativistic 
Thomson scattering theory Q 



77 max' 
E 9 
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We can see that the incident photon energy E p is boosted 
by a factor of approximately 4 7 2 after the backscatter- 
ing. Therefore, the Compton scattering of photons with 
relativistic electrons can be used to produce high energy 
photons, i.e., gamma-ray photons. 

Under a set of conditions 6*.; ss it and Of « 0, the uncer- 
tainties of the scattered photon energy E g due to the un- 
certainties of the variables (E ei E pi Of and 0i) in Eq. ([2|) 



can be estimated [lfj, . For example, the relative un- 
certainty of the scattered photon energy AE g /E g due to 
the uncertainty of the electron beam energy AE e /E e is 
given by taking the derivative of Eq. @ with respect to 
E e . i.e., 



E n 



2(1- 



2-f 2 E p /E e AE e AE e 



1 + A^Ep/E/ E, 



E, 



(7) 



Contributions to AE g /E g associated with other variables 
are summarized in Table |U 



B. Scattering cross section 



Lorentz invariant form 



The general problem concerning the collision is to find 
the probabilities of final states for a given initial state 
of the system, i.e., the scattering cross section. Using 
Quantum Electrodynamics (QED) theory, the Compton 
scattering cross section in the Lorentz invariant form has 
been calculated in [4j, ll2j, |l3j , and the result for unpolar- 
ized electrons scattering with polarized photons is given 
by 
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where r e is the classical electron radius; <pf is the az- 
imuthal angle of the scattered photon; £1,2,3 an d 2 3 
are Stokes parameters describing the incident and scat- 
tered photon polarizations in their respective coordinate 
systems; X and Y are the Lorentz invariant variables 
defined as follows 



X 



s — (mc)' 



Y 



(mc) 2 — u 



(9) 



(mc) 2 ' (mc) 2 
where s and u are the Mandelstam variables [3| given by 
8= (p + k) 2 , u = (p-k') 2 . (10) 
X and Y satisfy the inequalities [|[ 
X 



X + l 



< Y < X. 



(11) 



Since the scattering cross section of Eq. flSJ) is expressed 
in the Lorentz invariants, it can easily be expressed in 
terms of the collision parameters defined in any specific 
frame of reference. 



2. Polarization description in lab frame 

In the laboratory frame, three right-hand coordinate 
systems are used in Eq. flSJ) to describe the motion and 
polarization of the incident electron (x e ,y e , z e ), the inci- 
dent photon (x, y, z), and the scattered photon (x' , y', z') 
(Fig. [3]). The coordinate system (x e , y e , z e ) is fixed in the 
lab frame, and its z e -axis is along the incident direction 
of the electron, (x, y, z) and (x', y', z') are the local coor- 
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FIG. 3: Coordinate systems of Compton scattering of an elec- 
tron and a photon in a laboratory frame. (x e ,y e ,z e ) is the 
coordinate system for the incident electron (p) moving along 
the 2 e -axis direction. For the head-on collision, the incident 
photon (k) comes along the negative z e -axis, and the scattered 
photon (k 1 ) is moving along the direction given by the polar 
angle Of and azimuthal angle <f)f. The momentum vectors 
k and k' form the scattering plane, (x, y, z) is a right-hand 
coordinate system attached to the scattering plane. The z- 
axis is along the direction of k; 5-axis is perpendicular to the 
scatter plane; and y-axis is in the scattering plane, (x 1 , y , z') 
is another right-hand coordinate system attached to the scat- 
tering plane. The z'-axis is along the direction of k'; z'-axis 
is the same as the i-axis; and y'-axis lies in the scattering 
plane. 



dinate systems attached to the scattering plane formed 
by the momenta of the incident and scattered photons, 
k and k' . For (x, y, z), the x-axis is perpendicular to the 
scattering plane; the y- and z-axes are in the scatter- 
ing plane with the z-axis along the direction of k. For 
(x',y',z'), the i'-axis is the same as the i-axis for the 
incident photon, perpendicular to the scattering plane; 
and the z'-axis is along the direction of k! . 

The Stokes parameters ^\ 3 of the incident and scat- 
tered photons in Eq. ([5]) are defined in their local coordi- 
nate systems, respectively. The parameter £3 describes 
the linear polarization of the photon along the x^- or 
2/M-axis; the parameter £^ describes the linear polariza- 
tion along the direction at ±45° angles relative to the 
5;W-axis; and the parameter Q represents the degree of 
circular polarization of the photon. 

The polarization of the photon is always defined in its 
local coordinate system with its momentum being one of 



the axes. For Compton scattering described by Eq. ([8]), 
these local coordinate systems (x, y, z) and (x' , y',z') are 
different for different scattering planes. However, for the 
cases that the photons and electrons collide nearly head 
on to produce high-energy photons with small scatter- 
ing angles, it becomes possible to conveniently express in 
an approximate manner the polarization of the incident 
and scattered photons using a fixed coordinate system, 
for example, the lab-frame electron coordinate system 

(*£e; ye: 

Let us consider the incident photon with its z-axis ap- 
proximately parallel to the negative z e -axis. The Stokes 
parameter of the incident photon can be related to the 
degrees of polarization defined in the fixed electron coor- 
dinate system through the following equations [3] , 

6 w P t sin(2r-20 / )), 

£2 « Pc, 

£3 « -P t cos(2r-20 / )) 1 (12) 

where P t and P c arc the degree of linear and circular po- 
larizations of the incident photon defined in the coordi- 
nate system (x e ,y e , z e ), respectively; r is the azimuthal 
angle of the linear polarization P t with respect to the 
x^-axis ; and <pf is the azimuthal angle of the scattering 
plane. 

For Compton scattering involving an ultra-relativistic 
electron, scattered photons are concentrated in a small 
scattering angle (Of < 1/7). For these high-energy pho- 
tons with small scattering angles, their z'-axes are ap- 
proximately parallel to the z e -axis. Neglecting the polar 
angle (i.e. Of <C 1), the Stokes parameters of the scat- 
tered photon can be expressed approximately using a set 
of Stokes parameters defined in the fixed electron coor- 
dinate system as (l4| . 

£ ps -gcos20/+&sin2^/, 

£ 3 ps -£[ sin 20/ cos 20/, (13) 

where 2 3 are the Stokes parameters defined in the co- 
ordinate system [x e , y e , z e ). 

C. Spatial and energy distributions of scattered 
photons 

Based upon Eqs. ([SJ, (|12p and (fj"3")) . we can calcu- 
late the spatial and energy distributions of a gamma-ray 
beam produced by Compton scattering of a monoener- 
getic electron and laser beams with zero transverse beam 
sizes, i.e., the particle-particle scattering. 

Let us consider Compton scattering of an unpolarized 
electron and a polarized laser photon without regard to 
their polarizations after the scattering. The differential 
cross section is obtained by setting 2 3 to zero in Eq. © 
and multiplying the result by a factor of two for the sum- 
mation over the polarizations of the scattered photons [H . 
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Thus, the differential cross section is given by [ll| 
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ir'i 




IP' 




fx 







1 1 
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Y 



(14) 



The total cross section can be obtained by integrating 
Eq. (fT4")l with respect to Y and 4>f, 



o 'tot 



X 



^ri- ill---— log(l + X 



X X 2 
1 



X 2(1 + X) 



(15) 



Note that the Stokes parameter £3 depends on </>/; how- 
ever, after integration over <pf the dependence vanishes. 
Neglecting the recoil effect (X <C 1), we can have 



a tot 



-(1-X) 



8irr 2 



(16) 



3 v ' 3 ' 
which is just the classical Thomson cross section. 

1. Spatial distribution 

For a head-on collision (Oi — n) in a laboratory frame, 
according to Eq. ([9]) the Lorentz invariant quantities X 
and Y are given by 

x = 2 7 £ p (l+/^ y = 2 7 £: g (1-/3 cos 



and 



dy = 2 



sin 0fd9f. 



(18) 



Substituting dF into Eq. (|14"|) . the angular differential 
cross section is given by 



da 8r 2 I 

^ = ^| [l + P t cos(2r-2^)] 



1 1 

X ~ Y 
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1 (X Y 

~4\Y + X 



E„ 



(19) 



where dfi = sin0/d0/d<^>/ and £3 has been expressed in 
terms of P t (Eq. $5$). 

From Eq. (fl"9]l . we can see that the differential cross 
section depends on the azimuthal angle <j)f of the scat- 
tered photon through the term P t cos(2t — 2(f) f). For a 
circularly polarized or unpolarized incident photon beam 
(Pt = 0), this dependency vanishes. Therefore, the dis- 
tribution of scattered photons is azimuthally symmetric. 
However, for a linearly polarized incident photon beam 
(Pt 7^ 0), the differential cross section is azimuthally 



modulated, and the gamma photon distribution is az- 
imuthally asymmetric. Figs.0]and[5]illustrate the spatial 
distributions of Compton gamma photons at a location 
60 meters downstream from the collision point for both 
circularly and linearly polarized incident photon beams. 
In these figures we can also see that the distribution of 
scattered photons peaks sharply along the direction of 
the incident electron beam. This demonstrates that the 
gamma-ray photons produced by Compton scattering of 
a relativistic electron beam and a laser beam are mostly 
scattered into the electron beam direction within a nar- 
row cone. 



2. Energy distribution 

For a head-on collision in the laboratory frame, it can 
be shown that 



Y = X 



f3E e — Eg 

f3E e — E p 



Thus, 



dY 



-X- 



dE n 



j3E e — E p 



(20) 



(21) 



Substituting dy in Eq. (fl4l and integrating the result 
with respect to the azimuth angle 4>f, we can obtain the 
energy distribution of scattered photons as follows 
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dE~ n 
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(22) 



The energy spectrum calculated using Eq. (j22j) is 
shown in Fig. [51 The spectrum has a high energy cutoff 
edge which is determined by the incident electron and 
photon energies according to Eq. ([5]). In Fig. [6j we can 
see the spectral intensity has a maximum value at the 
scattering angle Of = 0, and a minimum value around 
the scattering angle Of — 1/7. The ratio between them 
is about 2 when the recoil effect is negligible. This will 
be shown in the next section. 

Note that the energy spectrum shown in Fig. [S] is for 
a Compton gamma-ray beam without collimation. How- 
ever, if the gamma-ray beam is collimated by a round 
aperture with a radius of R and distance L from the col- 
lision point, the energy spectrum will have a low energy 
cutoff edge, and its value can be calculated using Eq. ^ 
with Of = R/L. 



3. Observations for a small recoil effect 

For a small recoil effect (X < 1), we can approximate 
Eqs. (flU]) and (f2"2")l to draw several useful conclusions. 
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FIG. 4: (Color) The computed spatial distribution of Compton gamma-ray photons produced by a head-on collision of a 
circularly polarized 800 nm laser beam with an unpolarized 500 MeV electron beam. The distribution is calculated for a 
location 60 meters downstream from the collision point. The left plot is a 3-dimensional intensity distribution, and the right 
plot is the contour plot of the gamma-beam intensity distribution. 




FIG. 5: (Color) The computed spatial distribution of Compton gamma-ray photons produced by a head-on collision of a linearly 
polarized 800 nm laser beam with an unpolarized 500 MeV electron beam. The polarization of the incident photon beam is 
along the horizontal direction. The distribution is calculated for a location 60 meters downstream from the collision point. 
The left plot is a 3-dimensional intensity distribution, and the right plot is the contour plot of the gamma-beam intensity 
distribution. 



For convenience, we first define 

f(y)-(---) 2 + --- + -(- + -)- (23) 

JK ' \X YJ X Y 4 \Y X J V ' 
Using the inequality Eq. (fTTj) . it can be found that 

1 , . 2 + X , , 



approximately (with a negligible recoil effect, X <C 1), 

\<f{Y)<\. (25) 

Thus, the maximum and minimum spectral flux of the 
Compton gamma-ray beam are given by 

do_ _ 8nr 2 e 2 + X 
dE g nax ~ X((3E e - E p ) 4 ' { ° J 
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2 3 4 5 
Gamma-ray energy (MeV) 

FIG. 6: The computed energy distribution of Compton 
gamma-ray photons produced by a head-on collision of a 
800 nm laser beam with a 500 MeV electron beam. The scaled 
scattering angle jOf with the electron Lorentz factor versus 
the gamma-ray photon energy is also shown in the plot. The 
solid line represents the energy distribution of the gamma-ray 
photons, and the dash line represents the relation between the 
scaled scattering angle and photon energy. 



For a circularly polarized or unpolarized incident pho- 
ton beam, according to Eq. (fll?|) . it can also be calculated 
that the angular intensity of scattered gamma-ray pho- 
tons at the scattering angle Of = I/7 is about 1/8 of the 
maximum intensity at the scattering angle Of = 0, i.e., 



(d<7/dfi) e/=0 ~ 8' 



(32) 



In addition, integrating Eq. (fT4")l over the entire solid 
angle of the cone with a half-opening angle of I/7, i.e., 
integrating Y over the range of X(l — X/2) ^ Y ^ X(l — 
X) and 4>f over the range from to 2n, we can have 
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, do 8ir 
{ dE~/ min = XljE e 



E p )4(l + X) 



The ratio between them is 

{da/dE g )mg X 
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(28) 



which is shown in Fig. [6] 
When Of = 0, we can have 



£ 9 « 4 7 2 £ p , Y*X{1-X). (29) 

Substituting Y in Eq. J23J), we have f(Y) sa 1/2. Thus, 
the spectral flux has a maximum value around the scat- 
tering angle Of = 0. When Of — I/7, we can have 



E a 



X(l-f). 



(30) 



Substituting Y into Eq. d23J, we have f(Y) w 1/4. 
Therefore, the spectral flux has a minimum value around 
the scattering angle Of = I/7. These results are illus- 
trated in Fig. |51 

Expressed in terms of the total scattering cross sec- 
tion of Eq. (fT6|) . the fraction of scattered photons in the 
energy range [E g nax — AE g nax , E' g nax ] can be found ap- 
proximately as 



3(2 + X) AE ; 



■max 
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/\ jprnax 
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(31) 



This is a simple formula which can be used to estimate 
the portion of the total gamma-ray flux with a desirable 
energy spread AE" 1 ^ after collimation. 



Comparing Eq. (|3"31 to the total cross section of Eq. p^|) , 
we can conclude that about half of the total gamma-ray 
photons are scattered into the I/7 cone. This can be 
explained by considering the Compton scattering in the 
electron rest frame. In this frame, the Compton scatter- 
ing process is just like "dipole" radiation: the gamma- 
ray photons are scattered in all the directions, a half 
of the gamma photons is scattered into the forward di- 
rection, and the other half into the backward direction. 
When transformed to the laboratory frame, the gamma- 
ray photon scattered into the forward direction in the 
rest frame will be concentrated in the 1/7 cone in the 
laboratory frame. 



D. Polarization of scattered photons 



For polarized photons scattering with unpolarized elec- 
trons without regard to the final electron polarization, 
the cross section is given by Eq. ©. Substituting £1,2,3 
and £1 2 3 using Eqs. (IT2l and (fH?]) . and assuming the lin- 
ear polarization of the incident photon beam is along the 
x e -axis, i.e., r = 0, we can get 
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It should be noted that the Stokes parameters £j 2 3 
describe the polarization of the scattered photon selected 
by a detector, not the polarization of the photon itself 0. 
In order to distinguish them from the detected Stokes 
parameters £123; we denote the Stokes parameters of 
the scattered photon itself by £123- According to the 

/ 



rules presented in section 65 of [4j, £[ 2 3 are given by 
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(36) 



Integrating Eq. (|34|) over the azimuthal angle 0/ gives 
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Therefore, the averaged Stokes parameters of the scat- 
tered photons over the angle 0/ are given by (q) — 
($j)/($o), which depend on the incident photon polar- 
ization and variables X and Y. 

For example, for 100% horizontally polarized (P t = 
1,P C = 0, t = 0) incident photons scattering with unpo- 
larized electrons, the average Stokes parameters of the 
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FIG. 7: The average Stokes parameter (£g) of Compton 
gamma-ray photons produced by a 100% horizontally polar- 
ized (Pt = 1,P C — 0, r = 0) 800 nm laser photons head-on 
colliding with an unpolarized 500 MeV electrons. 



scattered photons are given by 

<*l) 



(*0> 
(*3> 



0, 



(£2 



( $ o) 4(4-A)' + 4-A + $ + $ 



(39) 



Clearly, the scattered photons retain the polarization of 
the incident photons. (£3) as a function of the scattered 
photon energy is shown in Fig. [7] for 800 nm laser pho- 
tons head-on colliding with 500 MeV electrons. It can be 
seen that the average Stokes parameter (£3} of scattered 
gamma-ray photons is almost equal to 1 around the max- 
imum scattered photon energy as in this case the recoil 
effect is negligible. It means the scattered gamma-ray 
photons with the maximum energy are almost 100% hor- 
izontally polarized. 



III. BEAM-BEAM SCATTERING 

In the previous section we discussed the spatial and 
spectral distributions of a gamma-ray beam produced 
by Compton scattering of monoenergetic electron and 
laser beams with zero transverse beam sizes, i.e., particle- 
particle scattering. However, in the reality, the incoming 
electron and laser beams have finite spatial and energy 
distributions, which will change the distributions of the 
scattered gamma-ray beam. Therefore, there remains 
a need to understand the characteristics of a Compton 
gamma-ray beam produced by scattering of a laser beam 
and an electron beam with specific spatial and energy 
distributions, i.e., the beam-beam scattering. 

In this section, we discuss the beam-beam Compton 
scattering process. First, we derive a simple formula to 
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FIG. 8: Compton scattering of a pulsed electron beam and a 
pulsed laser beam in the laboratory frame. Two coordinate 
systems are defined to describe electron and laser beams: the 
first coordinate system (a;, y, z) is the electron-beam coordi- 
nate system in which the electron beam is moving along the 
z-axis direction; the (xi,yi,zi) system is the laser-beam co- 
ordinate system in which the laser beam propagates in the 
negative z;-axis direction. The coordinate systems (x,y,z) 
and (xi,yi,zi) share the same origin. 



calculate the total flux of the Compton gamma-ray beam. 
Then, we present two methods, a semi-analytical calcu- 
lation and a Monte Carlo simulation, to study the spa- 
tial and spectral distributions of the gamma-ray beam. 
Based upon these methods, two computing codes, a nu- 
merical integration code and a Monte Carlo simulation 
code have been developed. These two codes have been 



benchmarked against the experimental results at High In- 
tensity Gamma-ray Source (HI7S) facility at Duke Uni- 
versity. 



A. Geometry of beam-beam scattering 

Figure [8] shows Compton scattering of a pulsed electron 
beam and a pulsed laser beam in a laboratory frame. Two 
coordinate systems are used: (x, y, z) for the electron- 
beam moving along the z-direction; the (xi,yi,zi) for 
the laser-beam propagating in the negative z;-direction. 
These two coordinate systems share a common origin. 
The time t = is chosen for the instant when the cen- 
ters of the electron beam and laser pulse arrive at the 
origin. The definition of these two coordinate systems 
allows the study of the Compton scattering process with 
an arbitrary collision angle, i.e, the angle between z-axis 
and negative z;-axis. For a head-on collision, the colli- 
sion angle equals 7r. In this case, the electron and laser 
coordinate systems coincide. 

In these coordinate systems, the electron and laser 
beams with Gaussian distributions in their phase spaces 
can be described by their respective intensity functions 
as follows M 



f e (x,y,z,x',y',p,t)-- 
f P (xi,yi,zi,k,t)-- 



1 



(2Tr) 3 e x e y a p a, 
1 



-exp 



47T 2 er;cr fc cr£ 



■exp 



j x x 2 + 2a x xx' + /3 x x' 2 j y y 2 + 2a y yy' + j3 y y' 2 (p - p Q ) 2 (z - ctf 
j + yf (zi + ct) 2 (k - k Q f 



2e % 



2a 2 



2(7? 



2a, 2 



2a 2 



2°2 



4-7T 



1 + ^7 



(40) 



p is the momentum of an electron, and po is the centroid 
momentum of the electron beam; x 1 and y' are the an- 
gular divergences of the electron beam in the x- and y- 
directions, respectively; a XtV , p x ,y and J x , y are Twiss pa- 
rameters of the electron beam; a p , a z and s x>y are the 
electron beam momentum spread, RMS bunch length, 
and transverse emittance, respectively; k and A are the 
wavenumber and wavelength of a laser photon, and ko is 
the centroid wavenumber of the laser beam; /3o,Cfc and 
a 1 are the Rayleigh range, the RMS energy spread and 
bunch length of the laser beam. Note that the waist of 
the laser beam is assumed to be at the origin of both 
coordinate systems. 



B. Total flux 

The number of collisions occurring during a time dt 
and inside a phase space volume d 3 p d 3 k dV is given 



by 1 

dN(r,p,k,t) = a tot (p,k)c(l- P-k/\k\)n e (r,p,t) 

xn p (f,k,t)d 3 p d 3 k dVdt, (41) 

where a tot (j>> k) is the total Compton scattering cross 
section which is determined by the momenta of the 
incident electron and laser photon, p and hk; /3 = 
v e /c is the relative velocity of the incident electron; 
n e (r,p,i) = N e f e (f,p,t) and n p (r,k,t) = N p f p (f,k,t), 
where / e (r, p, t) and f p (r,k,t) are the phase space in- 
tensity functions of electron beam and laser pulse, and 
iV e and N p are the total numbers of electrons and laser 
photons in their respective pulses. 

To calculate the total number of scattered gamma-ray 
photons produced by collision, Eq. (l4"Tj) needs to be inte- 
grated for the entire phase space and the collision time, 
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X A 



Ntot = J dN(r,p,k,t) 

= N e N p J a tot {p,k)c(l - pcosdi) 

x/ e (r*,p J t)/ p (f ? ,fc ) i)d 3 p d 3 fc dVAt. (42) 

where 9i is the collision angle between the incident elec- 
tron and laser photon. Assuming collisions occur at the 
waists of both beams (a x — a y — 0, a w — ^ X/3q / (4tt)) , 
the spatial and momentum phase space in the density 
functions can be separated, i.e, fei^P^t) = fe{^^)fe(p) 
and f p (f,k,t) = f P (r,t)f p (k). Since the cross section 
ctotip, k) only depends on p and k, we can have 

N tot = N e N p J C sc a tot (p, k)f e (p)f p (k)d 3 p d 3 fc, (43) 

where 

£ sc = c(l-£cos0i) J f e (r,t)f p (r,t)dVdt (44) 

is the single-collision luminosity defined as the number 
of scattering events produced per unit scattering cross 
section, which has dimensions of 1/ area [ll| . For a head- 
on collision (9i = n) of a relativistic electron (/? ss 1) and 
a photon, the single-collision luminosity can be simplified 
to 



Thus, Eq. (|43|) can be rewritten in a simple form 



Nt, 



N e N p C sc a t c 



(45) 



(46) 



where ot^t is the total Compton scattering cross sec- 
tion averaged over the momenta of the incident electrons 
and photons. Neglecting the energy spread of the elec- 
trons and photons, Otot can be approximated by Otot of 
Eq. (|15p . which can be further simplified to the classical 
Thomson cross section if the recoil effect is negligible. 

If the beam-beam collision rate is /o, the gamma-ray 
flux is given by 



dN tl 



dt 



N e N p C sc a to tfo- 



(47) 



C. Spatial and energy distributions: 
semi-analytical calculation 

To obtain the spatial and energy distributions of a 
Compton gamma-ray beam, the differential cross sec- 
tion should be used instead of the total cross section in 
Eq. (l4"2l . In addition, two constraints need to be imposed 
during the integration of Eq. (|4"2"T) [9, 10]. 



collimation plane 

detection 
point 




t> z 



FIG. 9: Geometric constraint for a scattered gamma-ray pho- 
ton. Fhe diagram only shows the projection of the constraint 
in the x-z plane. 



First, let us consider the geometric constraint, which 
assures the gamma-ray photon generated at the location 
r can reach the location fd shown in Fig. [9] In terms of 
the position vector, this constraint is given by 



\k'\ ~ \fa-A'' 



(48) 



where k' represents the momentum of the gamma-ray 
photon; r — (x, y, z) denotes the location of the collision; 
and fd = (xd, yd, Zd) denotes the location where the scat- 
tered gamma-ray photon is detected. Due to the finite 
spatial distribution and angular divergence of the elec- 
tron beam, a gamma-ray photon reaching the location 
fd can be scattered from an electron at different collision 
points with different angular divergences. 

The constraint of Eq. (|48l) projected in the x-z and y-z 
planes is given by 



Xd - x 



y 



Vd-y 
L 



(49) 



Here, X and 9 y are the projections of the scattering an- 
gle 9f in the x-z and y-z planes, i.e., 6 X — 8fCOS(j>f, 
By = 9fsm(j)f and 8^ — 9 X + 9y, where 9f and <f)f 
are the angles defined in the electron coordinate system 
(x e , y e ,Ze) in which the electron is incident along the Re- 
direction (Fig. [H]). x' and y' are the angular divergences of 
the incident electron, i.e., the angles between the electron 
momentum and z-axis. L is the distance between the col- 
lision point and the detection plane (or the collimation 
plane). Note that a far field detection (or collimation) 
has been assumed, i.e., L ^> |f] and L rj \?d\- 

The second constraint is the energy conservation. Due 
to the finite energy spread of the electron beam, the 
gamma-ray photon with an energy of E g can be produced 
by electrons with various energies and scattering angles. 
Mathematically, this constraint is given by 



S(E g -E g ), 



(50) 
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where 



E 9 = 



1 + i 2 9 f + A^Ep/mc 2 



(51) 



Imposing the geometric and energy constraints in 
Eq. (22), the spatial and energy distributions of a Comp- 
ton gamma-ray beam can be obtained by integrating all 
the individual scattering events, i.e., 

dN(E g ,x d ,y d ) f da - 

dn d dE g * N e N p J —6(E g -E g )c(l + P) 

xf e (x,y,z,x',y',p,t) 

x f p (x, y,z,k, t)dx'dy'dpdkdVdt, (52) 

where dfl d — dx d dy d /L 2 , and der/dil is the differential 
Compton scattering cross section. Note that a head-on 



collision between electron and laser beams has been as- 
sumed, and the density function f e (r,p,t) has been re- 
placed with f e (x,y, z,x' ,y' ,p,t) of Eq. (|40|) under the 
approximation p z ~ p for a relativistic electron beam. 
In addition, the integration J ■■■ f p (r, k, t)d 3 k is replaced 
with J • • • f p (x,y,z,k,t)dk, where f p (x,y, z,k,t) is de- 
fined in Eq. (|40[) . Integrations over dfc^ and dk y have 
been carried out since the differential cross section has 
a very weak dependency on k x and k y for a relativistic 
electron beam. 

Assuming head-on collisions for each individual scat- 
tering event (9i = n and da/dtt is given by Eq. (TIT))) ), 
neglecting the angular divergences of the laser beam and 
replacing x' and y' with 9 X and 9 y , we can integrate 
Eq. (|5"2"|) over dV, dt and dp to yield the following result 
(see Appendix E|, 



dN(E g ,x d ,y d ) r 2 e L 2 N e N p f°° fV*^. 



dE g dx d dy d 



where 



47r 3 Sc/3 o- 7 cr fc J Q J-^/ iEp /E g J-e xmax ^fU^oexOey 1 + 2jE p /mc 2 



4~f 2 E p E g {l +1 2 9 2 ) 



x exp 



E g {l+ 1 2 9 2 ) 

{6 X - x d /Lf 
2al 



C = i + K-^) 2 



2k/3 x e x 

2 k f3y £ y 



1% 



^ 2 E p 



2cos 2 (r - (f>f) 



(l+7 2 ^ 2 ) 2 



y d /L) 2 ( 7 - 70 ) 2 (fc-fco) 2 



2a 2 



2o* 



d6 x d6 y dk, 



(53) 



, Cx 



2kj3 x s x 

1 H a , &0x 



2 k /3y £ y 



PxZ' 

I £ y£,y 



ftyCy 



= J8i 



n w xmax 



AE p /E g - 9 2 er 7 = — |, 
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(54) 



and ffB e is the RMS energy spread of the electron beam. 

In a storage ring, the vertical emittance of the electron beam is typically much smaller than the horizontal emittance. 
For a Compton scattering occurring at a location with similar horizontal and vertical beta functions {fi x ~ /3 y ), the 
vertical divergence of the electron beam can be neglected. In addition, the photon energy spread of a laser beam is 
small, and its impact can also be neglected in many practical cases. Under these circumstances, the cross section 
term in Eq. (|53[) has a weak dependence on 9 y (ss y d /L) and k (w ko). With the assumption of an unpolarized or 
circularly polarized laser beam, Eq. (|53")) can be simplified further after integrating 9 y and k: 



dN(E g ,x d ,y d ) 
dE g dx d dy d 



r 2 L 2 N P N n 



-yV(j x J- 



1-9) 



j 2 9 2 ) 2 X 6XP 



1 



1 + 2-fE p /mc 2 I 4 



4 7 2 £ P , E g {l + 1 2 9 2 ) 



E g (l + 1 2 9 2 ) 



x d /L) 2 ( 7 - 7o ) 5 



2 < 



2a 2 



d9 x , 



4 7 2 £ p 



(55) 



where 9 , 



= y/4E p /E g - {y d /Lf 



The integrations with respect to k, 9 y and 9 X in 
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Eq. (153)) or 6 X in Eq. (|55[) must be carried out numeri- 
cally. For this purpose, a numerical integration Compton 
scattering code (CCSC) in the C++ computing language 
has been developed to evaluate the integrals of Eqs. (|53")l 
and (|55i 

With the detailed spatial and energy dis- 
tributions of the Compton gamma-ray beam 
dN(E g ,x d ,y d )/{dE g dx d dy d ) 1 the energy spectrum 
of the gamma-ray beam collimated by a round aperture 
with a radius of R can be easily obtained by integrating 
dN(Eg,x d ,y d )/(dEgdx d dy d ) over the variables x d and 
y d for the entire opening aperture, i.e., \J x d + y\ ^ R 2 . 

The transverse misalignment effect of the collimator 
on the gamma-ray beam distributions can be introduced 
by replacing x d and y d with x d + Ax and y d + Ay in 
Eq. (l53l) or Eq. (|55|) . where A:r and Ay are the collima- 
tor offset errors in the horizontal and vertical directions, 
respectively. 



D. Spatial and energy distributions: Monte Carlo 
simulation 



In the previous section, we have derived an analyti- 
cal formula to study the spatial and energy distributions 
of a Compton gamma-ray beam. However, to simplify 
the calculation several approximations have been made: 
head-on collisions for each individual scattering event, a 
negligible angular divergence of the laser beam, and far 
field collimation. 

A completely different approach to study Compton 
scattering process is to use a Monte Carlo simulation. 
With this numerical technique, effects that cannot be 
easily included in an analytical method can be properly 
accounted for. For example, using a Monte Carlo simula- 
tion we can study the scattering process for an arbitrary 
collision angle. With this motivation, we developed a 
Monte Carlo Compton scattering code. In following, the 
algorithm of this code is presented. 



1. Simulation setup 

At the beginning of the collision, both the electron 
and laser pulses are located some distance away from 
the origin (Fig. [5]), and two pulse centers arrive at the 
origin at the same time (t = 0). The collision duration 
is divided into a number of time steps, and the time step 
number represents the time in the simulation. 

Due to a large number of electrons in the bunch, it 
is not practical to track each electron in the simulation. 
Therefore, the electron bunch is divided into a number 
of macro particles (for example, 10 6 ) which are tracked 
in the simulation. 

The phase space coordinates of each macro particle 
are sampled at time t = 0. For an electron beam with 
Gaussian distributions in phase space, the coordinates 



are sampled according to the electron beam Twiss pa- 
rameters as follows LL6L LL3 



x(0) 
x'(0) 
2/(0) 

y'(o) 

z(0) 
E„ 



\f2ui£ x p x cos <t>\, 

- y / 2u 1 £ x /f3 x (a x cos 0i + sin 0! ) , 

y/2u 2 Sy(3y cos </> 2 , 

- ^2u2£ v l (3 y {a v cos 02 + sin 4> 2 ) , 

E (l + cr_E c r 2 ), 



(56) 



where u\ t 2 are random numbers generated using an ex- 
ponential distribution with a unit mean parameter (i.e., 
e _tl1 ' 2 ), ri t 2 are random numbers generated according 
to a Gaussian distribution with a zero mean and unit 
standard deviation, and 4>i,2 & re uniformly distributed 
random numbers between and 2-7T. The coordinates of 
macro particles at any other time (t ^ 0) can then be ob- 
tained by transforming the coordinates given by Eq. (|5rJj) . 

The Compton scattering is simulated according to the 
local intensity and momentum of the laser beam at the 
collision point. The intensity of the laser beam at the 
collision point (x,y,z) in the electron-beam coordinate 
system can be calculated according to Eq. (t40|) using 
the laser-beam coordinates (xi,yi,zi) transformed from 
(x, y, z). The momentum direction k of the photon at the 
collision point (x, y, z) can be calculated from the point 
of view of electromagnetic wave of the photon beam. For 
a Gaussian laser beam, its propagation phase il)(xuyi, zj) 
in the laser-beam coordinate system is given by 



ip(xi,yi, z{) = -ikizi - ikiZi 



Vl 



(57) 



the wavevector (the momentum of photon ki) is given by 
h = yip(xi,yi,zi). Thus, 



l 



yji + c\ + C 



(cixi + c 2 yi + zi), (58) 



where 



Cl 



XlZi 
?2 , ^2 



f'2 



(59) 



The unit vector fc; expressed in the electron-beam coor- 
dinate system gives the momentum direction of the laser 
photon in this coordinate system. 



2. Simulation procedures 

At each time step, the Compton scattering process is 
simulated for each macro particle. The simulation pro- 
ceeds in two stages. In the first stage, the scattering 
probability is calculated using the local intensity and mo- 
mentum of the laser beam. According to this probability, 
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the scattering event is sampled. If the scattering hap- 
pens, a gamma-ray photon will be generated, and the 
simulation proceeds to the next stage. In the second 
stage, the energy and scattering angles (including the 
polar and azimuthal angles) of the gamma-ray photon 
are sampled according to the differential Compton scat- 
tering cross section. The detailed simulation procedures 
for these two stages are presented as follows. 



3. First stage: scattering event 

Since the energy and scattering angles of the gamma- 
ray photon are not the concern at this stage, the total 
scattering cross section is used to calculate the scatter- 
ing probability. According to Eq. (|3Tj). the scattering 

probability P(r,p, k, t) in the time step At for the macro 
particle at the collision point (x, y, z) is given by 

P(r,p,k,t) = a to t(P,k)c(l - P ■ k/\k\)n p (x,y,z,k,t)At, 

(60) 

where n p (x, y, z, k, t) and k are the local density and 
wavevector of the photon beam, respectively; cr tot (p, fc) 
is the total scattering cross section given by Eq. (Tl"5|) . 

According to the probability P(r,p,k,t), the scatter- 
ing event is sampled using the rejection method as fol- 
lows [l9l . [20l | : first, a random number r% is uniformly 

generated in the range from to 1; if < P(f,p,k,t), 
Compton scattering happens; otherwise the scattering 
does not happen, and the above sampling process is re- 
peated for the next macro particle. 



4- Second stage: scattered photon energy and direction 

When a Compton scattering event happens, a gamma- 
ray photon is generated. The simulation proceeds to the 
next stage to determine the energy and scattering angles 
of the gamma-ray photon. For convenience, the sampling 
probability for generating gamma-ray photon parameters 
is calculated in the electron-rest frame coordinate system 
(x' e ,y' e ,z' e ) in which the electron is at rest and the laser 
photon is propagated along the z^-axis direction. 

Since the momenta of macro particles and laser pho- 
tons have been expressed in the electron-beam coordi- 
nate system (x, y, z) in the lab frame, we need to trans- 
form the momenta to those defined in the electron-rest 
frame coordinate system (x' e , y' e ,z' e ). After transforma- 
tions, the sampling probability for generating the scat- 
tered gamma-ray photon energy and direction will be cal- 
culated as follows. 

In the electron-rest frame coordinate system 
(x' e ,y' e ,z' e ), according to Eq. J2]) the scattered pho- 
ton energy is given by 



1 



1 

~e7,. 



1 



+ — ? (l-cos0'), 



(61) 



where 0' is the scattering angle between the momenta of 
the scattered and incident photons; E' g and E' p are the 
energies of the scattered and incident photons, and E' is 
in the range of 



1 + 2E> Jm& 



< E'< 



(62) 



In the electron-rest frame coordinate system, we can 
simplify the Lorentz invariant quantities X and Y of 
Eq. (JT3J) to A = 2E' p /mc 2 and Y = 2E' g /mc 2 . As a 
result, the differential cross section is given by 
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( mc 2 




P t cos(2t 




/ mc 2 
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E ' g 





EL 



^ + ^r , (63) 



where r' is the azimuthal angle of the linear polariza- 
tion direction of the incident photon beam defined in the 
system (x' e ,y' e , z' e ), and <j)' is the azimuthal angle of the 
scattered photon. Note that the quantity Pt, the degree 
of linear polarization of the incident photon beam, is in- 
variant under Lorentz transformations. 

The scattered photon energy E' g and the azimuthal 
angle <fi' are sampled according to the differential cross 
section Eq. (|S3"f . Since Eq. (|6"3"j) depends on both E' and 
<j>' the composition and rejection sampling method [l9l . 
[20 ] is used to sample these two variables. To sample the 
scattered gamma-ray photon energy E' g , Eq. (f6"3"|) needs 
to be summed over the azimuthal angle <fi' and written 
as 



da 



9 mc ,„ 
d^ = 7n % ( + 
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(64) 



where 



f(E') = 
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2 + 2E' p /mc 
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(65) 



and < f(E' g ) < 1 for any E' g . Now, the scattered 
gamma-ray photon energy E' can be sampled according 
to f{E' g ) as follows: first, a uniform random number E' g 
is generated in the range given by Eq. (|62"T) . and r± in 
the range from to 1; if < f{E' g ), E' g is accepted, 
otherwise the above sampling process is repeated until 
E g is accepted. If E' g is accepted, the scattering angle 6' 
can be calculated using Eq. (|61|) . 

After the scattered gamma-ray photon energy E' g is 
determined, the azimuthal <j)' angle is sampled according 
to 



d 2 a da 

dE'^'dE^' 



(66) 
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Gamma-ray energy (MeV) 

FIG. 10: Compton gamma-ray beam energy spectra gener- 
ated using computer codes MCCMPT, CCSC and CAIN2.35. 
The stairs plot represents the spectrum simulated using the 
code MCCMPT, the dash line represents the spectrum cal- 
culated using the code CCSC, and the circles represent the 
one using the code CAIN2.35. The electron beam energy and 
RMS energy spread are 400 MeV and 0.2%, respectively. The 
electron beam horizontal emittance is 10 nm-rad, and the ver- 
tical emittance is neglected. The laser wavelength is 600 nm 
with negligible photon beam energy spread. The gamma-ray 
beam is collimated by an aperture with a radius of 12 mm 
located 60 meters downstream from the collision point. 




Gamma-ray energy (MeV) 

FIG. 11: Comparison between the measured and calculated 
energy spectra of a Compton gamma-ray beam. The solid 
line represents the calculated spectrum using the CCSC code, 
and the circles represent the measured gamma-beam energy 
distribution after removing the escape peaks and Compton 
plateau using a spectrum unfolding technique. The gamma- 
ray beam is produced by Compton scattering of a 466 MeV 
electron beam and a 790 nm laser beam at the HI7S facility. 
The RMS energy spread of the electron beam is 0.1%, and 
horizontal and vertical emittance are 7.8 and 1.0 nm-rad, re- 
spectively. The collimator with an aperture radius of 12.7 mm 
is placed 60 meters downstream from the collision point. 



After obtaining the gamma-ray photon energy E' g , and 
the angles 9' and </>' in the electron-rest frame coordi- 
nate system, we need to transform these parameters to 
those in the lab-frame coordinate system. In the mean- 
time, the momentum of the scattered electron is also 
computed. This electron can still interact with the laser 
photon in following time steps, which allows to correctly 
model the multiple scattering process between the elec- 
trons and laser photons. 



IV. BENCHMARK AND APPLICATIONS OF 
COMPTON SCATTERING CODES 

Based upon the algorithms discussed in Section III, we 
have developed two computer codes using the CH — h pro- 
gramming language: the numerical integration Compton 
scattering code CCSC and the Monte Carlo Compton 
scattering code MCCMPT. Below, we briefly discuss the 
benchmark and applications of these two codes. 



A. Energy Distribution 

Our Compton scattering computer codes MCCMPT 
and CCSC have been benchmarked against a well known 
beam-beam colliding code CAIN2.35 developed at KEK 



for International Linear Collider [l6(. The energy spec- 
tra of Compton gamma-ray beams generated using these 
three codes are shown in Fig. [lOj We can see that 
these three codes can produce very close results. In 
terms of computing time, the codes CCSC, MCCMPT 
and CAIN2.35 took about 10 minutes, 150 minutes and 
1200 minutes to generate these spectra using a single- 
core machine, respectively. Compared to the multi- 
purpose beam-beam colliding code CAIN2.35, the ded- 
icated Compton scattering codes CCSC and MCCMPT 
are much faster and easy to use. 

At the HI7S facility, the Compton gamma-ray beam 
is usually measured using a high-purity Germanium 
(HPGe) detector. Due to the non-ideal response of the 
detector, the measured spectrum has a structure of a 
full energy peak, a single and double escape peaks, and 
a Compton plateau. To unfold the measured energy 
spectrum, a novel end-to-end spectrum reconstruction 
method has been recently developed .21]. The compar- 



ison of the measured gamma spectrum and calculated 
spectrum using the CCSC code is shown in Fig. [TT] A 
very good agreement between them is observed. 

Using the Monte Carlo simulation code, we can study 
the Compton scattering process with an arbitrary colli- 
sion angle. The simulated spectra using MCCMPT are 
compared to those using CAIN2.35 in Fig. [12] Again, 
very good agreements are observed. It is clearly shown 
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FIG. 12: Compton gamma beam energy spectra for different 
collision angles 90°, 100°, 120°, 135° and 180°. These spec- 
tra are simulated using codes MCCMPT and CAIN2.35. The 
electron beam and laser beam parameters are the same as 
those in Fig. 1101 The solid lines represent the spectra simu- 
lated using the code MCCMPT, and the circles represent the 
spectrum simulated using the code CAIN2.35. 



that the gamma-ray beam produced by a head-on col- 
lision of an electron and a laser beams has the highest 
energy and flux. With a 90 degree collision angle, the 
maximum energy of the gamma-ray beam is only half of 
that for a head-on collision. 

The energy spread of a Compton gamma-ray beam is 
mainly determined by the degree of the collimation of 
the gamma beam, energy spread and angular divergence 
of the electron beam [2lj . The contributions of these 
parameters to the gamma-ray beam energy spread are 
summarized in Tabic |IJ In some literature [22l . |23| , a sim- 
ple quadratic sum of individual contributions was used 
to estimate the energy spread of the Compton scattering 
gamma-ray beam. Since the electron beam angular diver- 
gence and the gamma-beam collimation introduce non- 
Gaussian broadening effects on the gamma-beam spec- 
trum [2lJ , causing the spectrum to have a long energy tail 
(Figs. [10] and [UJ, the energy spread of the gamma-ray 
beam cannot be given simply by the quadrature sum of 
different broadening mechanisms. The realistic gamma- 
ray beam energy spread needs to be calculated from its 
energy spectrum, which can be done using either the nu- 
merical integration code CCSC, or a Monte Carlo simu- 
lation code, MCCMPT or CAIN2.35. 



and linear polarized incoming laser beams. For compar- 
ison, the measured spatial distributions of gamma-ray 
beams using the recentl y d eveloped gamma-ray imaging 
system at HI7S facility [24| are also shown in Fig. [ljS] It 
can be seen that for a circularly polarized incoming laser 
beam, the distribution is azimuthally symmetric; for a 
linearly polarized incoming laser beam, the gamma-ray 
beam distribution is asymmetric, and is "pinched" along 
the direction of the laser beam polarization. 

More applications of using CCSC and MCCMPT codes 
to study characteristics of Compton gamma-ray beams 
can be found in ED, H, [3 . 



SUMMARY 



To study characteristics of a gamma-ray beam pro- 
duced by Compton scattering of an electron beam and 
a laser beam, we have developed two algorithms: one 
based upon an analytical calculation and the other us- 
ing a Monte Carlo simulation. According to these algo- 
rithms, two computer codes, a numerical integration code 
(CCSC) and a Monte Carlo simulation code (MCCMPT), 
have been developed at Duke University. These codes 
have been extensively benchmarked against a beam-beam 
colliding code CAIN2.35 developed at KEK and measure- 
ment results at the High Intensity Gamma-ray Source 
(HI7S) facility at Duke University. Using these two 
codes, we are able to characterize Compton gamma-ray 
beams with various electron and laser beam parameters, 
arbitrary collision angles, and different gamma-beam col- 
limation conditions. 

In this work, the nonlinear Compton scattering pro- 
cess is not considered, and the polarization of the elec- 
tron beam is not taken into account. Although the po- 
larization of the gamma-ray beam has been calculated 
in Section II, this calculation is limited to the particle- 
particle scattering case. Further studies will be carried 
out to address these issues. 
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FIG. 13: (Color) Spatial distributions of Compton gamma-ray beams for different polarizations of the incoming laser beams. 
The gamma-ray beams are produced by Compton scattering of a 680 MeV electron beam and a 378 nm laser beam. The 
observation plane is about 27 meters downstream from the collision point. The upper plots are the simulated images using 
the MCCMPT code. The lower ones are the measured images. The left images are produced using a circularly polarized laser 
beam; the right images with the linearly polarized laser beam. 



Appendix A: Spatial and energy distributions of a Compton gamma-ray beam 



The spatial and energy distributions of a Compton gamma-ray beam produced by a head-on collision of an electron 
beam and a photon beam is given by 



dN(E g ,x d ,y d ) 
dn d dE g 



^ s i E g - E g )c(l + /3)n e (x, y, z, x', y',p, t)n p (x, y, z, k, t)dx'dy'dpdkdVdt, (Al) 



where dil d — dx d dy d / L 2 ; n e (x, y, z, x' , y' ,p, t) and n p (x, y, z, k, t) are the density functions of the electron and photon 
beams given by Eq. (Pf(J)) ; der/dfi is the differential cross section given by Eq. (fTI?)) . For head-on collisions, we can 
simplify the differential cross section to 



da_ 



8ri 



47 2 £ P , E g (l + j 2 9}) 
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7 2 ^ 



( E 9 
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Replacing x' and y' with 9 X and 8 y according to Eq. and neglecting the angular divergence of the laser beam 
at the collision point, we can integrate Eq. (|A1|) over dT^ and dt and obtain 



dN(E g ,x d ,y d ) 
dE g dx d dy d 



L 2 N e N p 
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1 da 
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(A3) 
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where 



L> PO 



PV\2 | 2kf3yEy _ 

l ' + A) ' Cy ~ 



2k/3 x s a 
_ So - ' 



/ ^yCy 



PyCy 



9f = \J0% + 9y, 9 X = Of cos0/, y = f sincj) f . 



(A4) 

Next, we need to integrate the electron beam momentum dp. It is convenient to change the momentum p to the 



scaled electron beam energy variable 7 = E e /(mc ), and rewrite the delta-function S(E g — E g ) as 



8{E g - E g ) = S{ 



4j 2 E p 



where 



is the root of 



1 + l 2 2 } + AqEpjmc 2 

_ 2E g E p /mc 2 ( 
4Ep — E g 2 \ 

E n = 



E g ) = -5(7 - 7) 



(1 +1 2 0) +4jE p /mc 2 ) 2 
% 1 Ep{l + 2 1 Ep/mc 2 ) 



AE p - E g 2 
4E 2 E g /(mc 2 ) 2 



A^Ep 



+ -f 2 2 + 4jE p /mc 2 



with the condition of < Of < W -gr 1 

Substituting Eqs. (|A2j) . (| A5|) into Eq. (|A3|) and integrating d-y, we can get 



dN{E gi x d ,Vd) 
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where 



= x lAEp/E g -0 2 y . 
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